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' Dirac leptogenesis (or Dirac neutrinogenesis), in which neutrinos are purely Dirac particles, 

is an interesting alternative to the standard leptogenesis scenario. In its supersymmetric version, 
' the modified form of the superpotential required for successful baryogenesis contributes new, 

, generically non-flavor-diagonal terms to the slepton and sneutrino mass matrices. In this work, 

we examine how current experimental bounds on flavor-changing effects in the lepton sector 
(and particularly the bound on /i ^ ej) constrain Dirac leptogenesis and we find that it is 
capable of succeeding with superpartner masses as low as ~ 100 GeV. For such light scalars 
Ph. and electroweakinos, upcoming experiments such as MEG are generically expected to observe 

^ ' signals of lepton flavor violation. 

> 

X ■ 1 Introduction 

When singlet fermions are not present in a given theory, Dirac leptogenesis ^ or Dirac 
neutrinogenesis, represents a very interesting alternative to the traditional leptogenesis scenario, 
which relies on the existence of heavy Major ana neutrinos. In Dirac leptogenesis, neutrinos are 
purely Dirac particles whose small but nonzero masses appear as ratios of dimensionful parameters 
in an effective field theory. It has been shown (S] that, in the context of split supersymmetry [IJ 

Dirac leptogenesis is a phenomenologically viable scenario capable of satisfying all relevant 
constraints from cosmology and neutrino physics as well as reproducing the observed baryon-to- 
photon ratio r] of the universe. Split supersymmetry is advantageous primarily for two reasons. 
The first of these is that very heavy gravitinos can easily evade the constraints that big bang 
nucleosynthesis (BBN) places on the post-infiationary reheating temperature; the second is that in 
Dirac leptogenesis the superpotential is extended by the addition of new terms carrying phases and 
nontrivial flavor structure. Dangerous contributions from these terms to flavor-changing processes 
are automatically safe in split supersymmetry, but must be treated with care when the scale of all 
the superpartners is near the electroweak scale. In the latter case (assuming that the constraints for 
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gravitino cosmology are satisfied) it is necessary to compute carefully the rates for flavor-violating 
processes in the lepton sector. 

The aim of this paper is to ascertain whether or not Dirac Leptogenesis is permitted by flavor 
violation constraints in supersymmetric theories with low-scale slepton masses. We will begin 
our inquiry by briefly reviewing the theoretical framework of Dirac leptogenesis and deriving the 
additional contributions to the slepton mass matrices to which its superpotential gives rise. We then 
turn to the calculation of rates for the processes /i ^ 67 and t ^ ji^. Finally, we apply the combined 
constraints from flavor violation and baryogenesis and assess the viability and predictability of Dirac 
leptogenesis when superpartner masses are at or around the weak scale. 



2 Dirac Leptogenesis and the Slepton Mass Matrices 



In Dirac leptogenesis, as in the traditional leptogenesis picture [S'iT, the conditions for successful 
baryogenesis are met by positing the existence of a heavy particle with CP-violating decays into 
leptons.^ The lepton number produced in these decays is then processed into a baryon number 
for the universe by electroweak sphaleron processes ^llj. In the supersymmetric version of Dirac 
leptogenesis, all of this is engineered via a specific set of modifications to the superpotential and the 
postulation of a few additional superfields. In addition to the matter fields of the MSSM, at least 
two massive vector-like pairs of chiral superfields $ and $ are required, as are three generations 
of right-handed neutrino superfield Na and an additional exotic superfield Xj whose function will 
be to acquire a scalar VEV. An additional symmetry, whose breaking will be responsible for late 
neutrino masses, is also imposed, and charges under it assigned so that the most general leptonic- 
sector superpotential is 

W 3 XiaNa^.Hu + hic^La^iX + M$,$i¥i + fiHuHd- (1) 

Here, A and h are (generally complex) coupling matrices and M$. are the masses of the heavy 
vector- like pairs, which, in order for successful baryogenesis to occur, are generically required to 
be 10^° GeV or larger.^ Once the heavy superfields and are integrated out, the resulting 
effective superpotential 

\ . h* 

^-S 3 ^T^xLpH^N^ + nHuHd (2) 

will yield a small but nonzero Dirac mass term for the neutrinos, provided that some mechanism is 

invoked to give the scalar component of x a VEV. Since A and h are complex, the decays of both the 

scalar and fermionic components of <I>i and <I>i, the lightest of the additional heavy superfields, will 

in general be CP-violating'^ and will result in two equal and opposite stores of lepton number, Lagg 

and Lr. The first of these is an aggregate lepton number stored in left-handed leptons, sleptons, 

^For recent variations on both supersymmetric and non-supersymmetric scenarios, see for example |nillU|. 
■^In general M#, A, and h are all complex matrices, but we can always choose to work in a basis where M$ is 
diagonal. 

^It can be shown that there is at least one nontrivial CP-violating phase in A and h. 
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and other fields in equilibrium with them, and is transformed into baryon number by sphaleron 
interactions. The second is stored only in right-handed neutrinos,^ which, being singlets under 
SU{2) X U{1)y, do not experience sphaleron effects and only couple to the other light matter fields 
through the effective neutrino Yukawa interaction given by © with x ~^ ix)- This interaction is 
suppressed by {x)/M^^ and therefore the time scale for the equilibration of left-handed and right- 
handed stores of lepton number can be quite late. If the effective neutrino Yukawas are sufficiently 
small, ^ the equilibration time scale will be much longer than the time scale at which sphaleron 
processes effectively shut off pp. When this is the case, a net baryon number for the universe will 
already have frozen in and will persist unaltered until present time. 

In addition to providing a mechanism for baryogenesis, Dirac leptogenesis holds some interesting 
implications for neutrino physics. Here, the squared neutrino masses are given by 

2 2 (or 3) 3 ^ 

\mu\lp= {v{x)sm(3^ ^ ^K^>^hh*ochjfi ■ (3) 

ij=l 7=1 *i *i 

In order to reproduce the observed mass splittings ^21 

Amii = (7.9^0:6) X lO-^eV^, \^ml,\ = {2.2t',l) x 10-3eV^ (4) 

where Am?^- = mf — m'j, and the angles in the Umns matrix, it is necessary to impose a few 
conditions on the coupling matrices A and h, which determine the matrix structure in Q. There 
are several ways of doing this, but we will focus on one particular, theoretically-motivated scenario 
called constrained hierarchical Dirac leptogenesis (CHDL) P], which produces a normal hierarchy 
among neutrino masses. 

In CHDL, the appropriate matrix structure for A, h, and the (diagonal) mass matrix M$ is 
obtained by requiring that A and h both be antisymmetric with C(l) entries and by relating M$ 
to the fermion Yukawas through a flavor symmetry. This structure ensures that the neutrino 
mass matrix one obtains after integrating out the heavy fields $j and $j corresponds to a normal 
neutrino hierarchy, provided that there is a hierarchy among the mass eigenstates and M^j. 
Of course this coupling structure also impacts baryogenesis, but it has been shown that in CHDL, 
the constraints from baryogenesis and neutrino physics can be satisfied simultaneously 

Baryogenesis and a small neutrino Dirac mass are not the only consequences of equation (^, 
however. Assuming that at some high scale M all the soft supersymmetry breaking terms are 
flavor diagonal and universal, the new terms of equation will induce potentially off-diagonal 
contributions to the slepton mass matrices m'j^j^ and iti\^. As in Major ana leptogenesis 

some of these contributions come from running the scalar soft masses from some high scale M 
(the Planck scale, the GUT scale, etc.) down to M$j. Assuming the Universality condition at the 
high scale and that all soft A-terms are equal to the relevant yukawa coupling multiplied by the 

*We assume that right-handed sneutrinos couple strongly enough to other fields in the theory (e.g. with left- 
handed sneutrinos in the event that x has a large F-term VEV) to equilibrate with them and thus contribute to Lagg- 
Otherwise they would contribute to Lr. 

^Recall that in this scenario small (Dirac) neutrino masses are obtained through small (effective) neutrino Yukawas. 
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universal soft supersymmetric mass m^, one can simply estimate the flavor violating corrections 
to the mass matrix by integrating the RGE equations iteratively jl31 114j . and generally obtain 
off-diagonal contributions to the slepton masses Sm'j^^ and 5m\j^\ 



If the F-term of x acquires a VEV {F^) (a corollary in most mechanisms via which its scalar 
component obtains a VEV), the effective theory superpotential yields yet another potentially off 
diagonal scalar mass term 

S<R = hlKp^^^vsmP (7) 

after electroweak symmetry breaking, which mixes left-handed and right-handed sneutrinos. 

In CHDL, where we have a specific flavor structure for the matrices A, h and M$, one has a 
specific prediction for flavor mixing among sleptons once the electroweak and hidden symmetries 
are broken. In order to examine the effect of these mixings, the full mass matrices for both the 
charged sleptons and sneutrinos must be taken into account. For simplicity, we will continue to 
assume that the leading soft breaking sector is flavor diagonal and universal with a common scalar 
mass TJis. The resulting additional contributions to the slepton mass squared matrices, given by 
equations ©, and ((7)), can thus be expressed in terms of the 3x3 submatrices Smj^^, Smj^j^, 
and Sm'j^j^ as 

H± = ^ - - - Srnl=\ u . (8) 

V ; y \ (Sml^Y ; 6ml, J 

The only contribution to the charged slepton mass squared matrix comes from 5m\j^, while the 
sneutrino mass squared matrix receives not only additional flavor mixings among left-handed and 
among right-handed sneutrinos, but also an effective A-term from 5m\^ which intermixes left- 
handed and right-handed sneutrinos. 

Once the matrices A and h have been fixed, up to an overall scaling parameter^ /, to give 
the correct leptonic mixing matrix UmnS: the remaining free parameters of the model are (x), 
M$. , and /. After the mass scale of the neutrinos is fixed, these parameters are not independent 
anymore, and are related by the constraint embodied in equation @. In CHDL, where the lightest 
®In CHDL A and h are antisymmetric, with 0(1) entries up to an overall scaling parameter f 

(9) 





where the values of 03 and 63 depend on the hierarchy between and Ali,^: when M^j/M^j = m^/mg, 63 can 

vary between 1.4 and 2.9, while can vary between 35 and 90; when M^^/M,s,-^ = 10, 63 can vary between 1.4 and 
2.9, while 03 can vary between 1.5 and 4.5 0)- 
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neutrino is very light relative to the other two, the masses of the two heavier neutrinos are 



ml^ « Aml = {7.9t'oi)xlO-^eV' (10) 
ml^ « Am^i = (2.21°:^) X lO'W. (11) 

Plugging in these values and using a-^ = 4.5 and 63 = 2.2 (the values most advantageous for 
baryogenesis with a chosen hierarchy M^^/M^-^ = 10) in the CHDL parametrization © of A and 
h, we find that the relation between the parameters (x), M^-, and /, obtained from the neutrino 
mass, is 

sin^ = 1.009 X 10"^^ (12) 

Using this relation, we can express the overall dependence of the slepton mass terms in (jSJ, @, 
and (O on the relevant mass scales in the theory. 



M. 

''RR J ^ 



5<R OC /2 OC (13) 



Note that the proportionality constants for the bottom two equations are not dimensionless: the 
ones associated with 5m\j^ and Smj^j^ each contain a factor of and have mass dimension [m]'^, 
while the one associated with 5m\^ contains a factor of {F^)v and has mass dimension [m]^. 



3 Flavor Violation 



The most stringent constraints on flavor violation in the lepton sector come from measurements 
of the branching ratios for flavor- violating decays and conversions of heavy leptons, such as ^ — > 67, 
r — > fi^, n eee and — > eA. The current experimental limits on the 2-body decay processes 
are ^ 

BR{fi^e-f) < 1.2x10"", (14) 
BR{t n-f) < 1.1 X 10~^ (15) 

In the near future, the MEG experiment ^H] is expected to improve the current experimental bound 
on ^ — > 67 by several orders of magnitude, to 0(10"^^ — 10~^^) or lower. Other related projects, 
such as PRIME J7] (sensitive to fiA — > eA conversion) , are expected to go online over the next few 
years. Projects have also been proposed ^S] that would lower the bound from r ^ /i7 to 0(10"^). 

The effective interaction leading to lepton flavor violating decays of the form £i ij'j, where 
£i and ij are charged leptons, can be written as 

I = ieme^ Ui{q - p)cj^pqf^ {A^Pl + A^Pr) Uj{p)e*{q), (16) 
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Figure 1: Feynman diagrams giving the two leading order contributions to the flavor-changing 
process ijj due to sneutrino (left diagram) and charged slepton (right diagram) mass mixings. 

where q and p are the momenta of the photon and the outgoing lepton ij respectively, and mi. is 
the outgoing lepton mass. The resulting decay rate is 

mj - c 7) = ^^d^^i' + 1^""!')- (17) 

The leading contributions to the amplitudes and appear at one loop level and are shown in 
figure n They involve both a sneutrino (and chargino) mass eigenstate and charged slepton (and 
neutralino) mass eigenstate running in the loop. These amplitudes have been computed in jl3j for 
a general MSSM scenario^ but for completeness we include them in Appendix A. 

In order to proceed further, it will be necessary to make a few assumptions concerning the 
supersymmetric model parameters. Those relevant to a discussion of lepton-sector flavor violation 
include the gaugino masses Mi and M2, the Higgs mass parameter //, the ratio of Higgs VEVs 
tan and the soft masses for the sleptons. In our analysis, we choose the values Mi = 160 GeV, 
M2 = 220 GeV, /i = 260 GeV, and tan /3 = 3, 10 and 30. As for the slepton soft masses we will, 
as previously mentioned, assume a common scale and examine what effect varying has on 
BR{^ 67) and BR{t fj,^). We will assume that the scale at which soft masses are universal is 
M = 2 X 10^6 GeV, though the results are not particularly sensitive to this choice. 

The results of our calculation are displayed in figure |2 In the left panel, we show exclusion 
contours in Af$^-(x) space for = 200 GeV. The areas below and to the right of the lower contour 
(the white region) are excluded by the experimental bounds given in dJ. We also include contours 
demarcating the region wherein baryogenesis can succeed, which have been updated from |3] to 
include the effects^ of processes second order in $1 and $1 (see Appendix El- Contours have also 
been computed for r — > ^7, but the constraints they imply for the theory are far weaker than 
those from fi — > 67. In the right panel, we show how varying the universal scalar mass affects the 
branching ratio for fi — > 67, which reaches a maximum when nis is around the weak scale. This 
is to be expected: when mg is much larger than the weak scale both the slepton mass-squared 

^In our case we need to add three right handed sneutrinos, but it is trivial to extend the result to include six 
sneutrino mass eigenstates instead of three. 

*We thank A. Strumia for pointing out to us the potential numerical importance of these thermalization processes 
involving gauge interactions |19l 1^ . 
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Figure 2: Exclusion plots combining constraints from both leptogenesis and flavor violation in the 
process — > 67. The left-hand panel shows exclusion contours in M^^-{x) space for a universal 
scalar soft mass = 200 GeV, with tan [3 = 10; the right hand panel shows the variation of the 
branching ratio BR{fj, 67) with respect to using tan /? = 3, 10 and 30. In both plots, we have 
taken Mi = 160 GeV, M2 = 220 GeV, and fi = 260 GeV. We have also assigned the x superfield 
an F-term VEV (F-^) = 10^ GeV. Such a large VEV results in large trilinear couplings between 
Higgs fields and sneutrinos and therefore induces potentially sizeable mixings between left-handed 
and right-handed sneutrinos after electroweak symmetry breaking. In each plot, the thick solid 
contours represent the current experimental bound on the branching fraction H14|): the dashed lines 
represent the expected future experimental bound of 1.5 x 10~^^ from MEG. The thin solid contour 
in the left-hand panel delimits the region allowed by leptogenesis constraints. 

eigenvalues and the flavor-violating terms scale like and the sneutrino and charged slepton 
mixing matrices asymptote to a constant value, while the branching ratio is still suppressed by the 
masses running in the loop; as decreases below the weak scale, Sm'j^j^ and Sm^^ go to zero and 
the slepton masses are dominated by flavor diagonal electroweak contributions. We also observe 
that, as in the SUSY see-saw case jl3l I14| . the flavor violation rate is quite sensitive to tan/3. 

In interpreting the results in figure |5J it is useful to note that in the regions of the plot near the 
exclusion contours (where the branching ratio for ^ 67 is quite low), flavor- violating effects will 
be small. It is therefore valid to use the mass-insertion approximation there and treat Sm'j^^ (left), 
5m\^, and 5m?^^ as small corrections to the slepton propagators. Since the sneutrino propagator 
can receive mass insertions from all three, we will focus our analysis on the partial amplitude 
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from the sneutrino- mediated process (the left diagram in figure^. In figure |31 we hst the leading 
contributions to this partial amplitude involving each of Sm'j^j^ (left diagram), 5m\j^, and dmj^j^. 
Since there is no coupling between leptons and right-handed sneutrinos, corrections from Sm'j^j^ 
and Sm'j^j^ only appear at second and third order in the mass insertion expansion, respectively. 
Therefore, if there is no substantial hierarchy among these three sets of mixing terms, the primary 
source of flavor-violation comes from mixings between left-handed sleptons. In the mass insertion 
approximation, the branching ratio for such processes can naively be estimated as 

BR{„ - ei) oc j;t—^ (18) 

and therefore contours of branching ratio in the — (x) plane correspond to contours of Smj^j^. 
According to equation 5m\j^ = ciM$^/(x), where ci is a dimensionless proportionality 

constant with dimension [m]^, so the exclusion contour associated with left-left mixing takes the 
form 

\nM^ = \n{x)+CLL, (19) 

where Cll = ~ lii(^"^iL/ci)- The oblique, upper exclusion contour in figure 12 which embodies 
the constraint from left-left mixing, is associated with this linear equation. 




Figure 3: Feynman diagrams for the leading-order process involving 5m\j^ (left diagram), and 
for the leading process involving Smj^j^ (right diagram), with sneutrinos running in the loop in 
the mass-insertion approximation. Note that any process involving Smj^j^ necessarily involves two 
mass insertions, and any one involving 5m?^^ (given by the diagram on the right with an additional 
Srnjjj^ insertion) necessarily involves three. 

While some hierarchy among the mass insertion terms is necessary for Sm?^^ and Smj^j^ to be 
relevant, there is no a priori reason why such a hierarchy should not exist. The Smj^j^ contribu- 
tion ((Tj) is proportional to {F^), which is essentially a free parameter. As mentioned above, (F-^) is 
not even relevant to Dirac leptogenesis per se, but appears as a common side-effect of mechanisms 
for assigning the x superfield a scalar VEV. Still, in many such mechanisms |^, the scale \/{Fj 
can potentially be quite large (10^ GeV or higher), and if this is the case, contribution from 5m\^ 
could potentially be as important as those from 5m\i^. Let us assume for a moment that this is 
the case and examine the constraints related to &ni\^ and 5m\^ together. In regions of figure [21 
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where (x) is small, we now have 



BRifi 67) oc " s = constant (20) 



along any exclusion contour. Equation (|13() tells us that, 

dml^ = C2-, (21) 
X 

where C2 has mass dimension [m]^. Therefore, when Smj^j^ <^ ^i^' associated contour is given 

by 

ln(x) = Clr, (22) 

where Clr = ln(c2/5m|^^). This equation explains the behavior of the contour in the left panel of 
figure[21when (x) is small and the oblique bound from left-left mixing abruptly ells into a horizontal 
line — the bound from left-right mixing. In general, the effect of increasing (F-^) is to push this latter 
bound upward, and for \/{F]J > 10^, the entirety of parameter space is excluded. The effect of 
right-right mixing is higher order still and only becomes relevant in regions of parameter space 
where iVf$^ is large and (x) is small — regions already excluded by left-left mixing. 

For the numerical analysis shown in Figure HI we have taken y/{F]J = 10'' GeV which in some 
parts of the (x) — M^-^ plane, induces substantial mixing between left-handed and right-handed 
sneutrinos. Sizeable mixing of this sort opens the intriguing possibility that the lightest sneutrino 
could be a legitimate dark matter candidate [221 HSl [2S HHl in Dirac leptogenesis, a possibility 
which would be interesting to investigate in the future. 

The primary message of figure |21 is that the combined constraints from leptogenesis and flavor 
violation do not exclude Dirac leptogenesis in theories with low-scale sfermion masses. In general, 
the latter set of constraints tend to rule out theories with exceptionally high masses for the decaying 
particles <I>i and $1 and exceptionally low values for (x)-^ Given current bounds on lepton flavor 
violation, there is still a sizeable region of parameter space within which the scenario succeeds, even 
when nis is as low as 100 GeV. It is of interest, however, that MEG and the next generation of 
lepton flavor violation experiments will be able to probe the vast majority of this region, and the 
data from these experiments will be crucial in determining the viability of Dirac leptogenesis with 
weak-scale superpartners. 

As was mentioned in section [21 some new symmetry must be posited in order to obtain late 
neutrino mass generation. In Dirac leptogenesis, neutrino masses are the result of the scalar 
component of the x superfield acquiring a VEV which breaks this new symmetry, producing a 
Goldstone boson or pseudo-Goldstone boson, depending on the nature of the symmetry. Constraints 
on such bosons arise from both BBN and cosmic microwave background (CMB) considerations |26j 

^Recall that the parameters (x), Mj>j and / enter in the effective neutrino Yukawa couphng given by|5| so when 
one of them is increased, we have to appropriately tune the other two to maintain the neutrino Yukawas at a fixed 
value. Such an adjustment, in turn, changes the different contributions from these parameters to off-diagonal entries 
in the slepton mass matrices. 
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as well as from the detection of abnormalities in the neutrino flux associated with supernova 
events ^7\, and they can become problematic (depending on the mass of the Goldstone boson) when 
the symmetry-breaking VEV is less than around 1 GeV. The value of (x) required by leptogenesis 
constraints is around 10 GeV, and thus the cosmological complications associated with breaking 
the additional symmetry necessary for Dirac leptogenesis do not pose any problem for the theory. 



4 Conclusion 



It has already been shown that Dirac leptogenesis stands as a phenomenologically viable 
alternative to the standard leptogenesis picture. In theories with comparatively light scalars, how- 
ever, flavor violation becomes a serious concern commonly solved by the Universality assumption. 
The Dirac leptogenesis superpotential gives rise to new mass terms for sneutrinos and charged 
sleptons which induce mixings between flavor eigenstates after the breaking of both the electroweak 
symmetry and the new symmetry responsible for late neutrino mass generation. Experimental limits 
on flavor violation in heavy lepton decays such as ^ — > 67 significantly constrain any theory which 
permits slepton flavor mixing, and in Dirac leptogenesis these constraints translate into bounds on 
the theory parameters (the mass of the heavy decaying particle) and (x) (the VEV of the 

exotic scalar field). Baryogenesis requirements also place significant constraints on both of these 
parameters, and thus the question as to whether leptogenesis can be made to work at all when 
scalar masses are light is a highly nontrivial one. 

In this work, we have shown that even when the masses of supersymmetric scalars are small, 
substantial regions of parameter space exist for which Dirac leptogenesis succeeds in producing a 
realistic baryon asymmetry for the universe while respecting current bounds on flavor violation. 
This is true even when the masses of supersymmetric particles are as low as ~ 100 GeV. Interestingly 
enough, for such light flelds, which presumably can be discovered at the LHC, it is generically 
predicted that experiments such as MEG should have enough sensitivity to observe flavor changing 
effects if Dirac leptogenesis is in fact responsible for the baryon asymmetry of the universe. 
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A Effective Couplings 



For completeness, we list here the results used in our analysis for lepton flavor violating 
processes. The amplitudes and in equation ((T7)) were computed in and, with a trivial 
extension to include three right-handed sneutrinos, are given by 



j^L _ ^(c),L _|_ ^(n),L 



and 



R 



(23) 



where the individual amplitudes A^'^^'^, A^")'^, A^'^'>'^, and A^")'^ are 



4 6 

— w — 

9'7r2 /-^ /-^ m? 



32vr2 ^ m? 

A=i x=i ix 



NIaxN^ax: 



6(1 - XAxY 
X (1 - 6xAx + 3x^js: + 2x\x - ^x\x InxAx) 



'i V 



"li,- (1 



1 ^ ^ 1 

A=1X=1 



2; AX; 



r(l - X^^ + 2XAX InxAx) 



(24) 



1 



6(1 - XAx] 



X (2 + 3xAx - ^x\x + Xax + 6xAX In xax) 

■3 + 4xAx - x\x - 2 In xax) 



R* 



Xa 



'lAX'-^-jAX 

mij (1 

A (n,c)L\ 
A'^ ' \l^r. 



XAX] 



(25) 
(26) 



Here, the indices A and X respectively label the gaugino (chargino or neutralino) and slepton 
(sneutrino or charged slepton) mass eigenstates, xax = "^xa/'^^x' ^I'ax i-^iAx) denote the 
effective couplings of charged lepton i to chargino (neutralino) A and sneutrino (charged slepton) 
X. The flavor mixing terms in (??) enter into the overall rate H17() through C^^^ and N^J^, 
which contain elements of the matrices Uu and that diagonalize the mass-squared matrices 
for sneutrinos and charged sleptons, respectively. The slepton masses also enter into the partial 
amplitudes (|24I26|) . 

The effective couplings N^^J^ and C^'J^ are 



R 

iAX 
L 

AX 



-SO 

92 f_ 



mw cos P 



{UN)A3Ux,i+3 I , 



riR 
'-'iAX 

'-'iAX 



y/2 ymw COS fi 
-92{Or)aiUxa, and 



{UN)A3Ui^, + 2{Un)ai ta^nOwU^x, 



i+3 ' 



92 



V2m\Y cos P 



{Ol)a2UI, 



(27) 



in terms of the chargino mixing matrices {Ofi)A.i and {OL)A,i the neutralino mixing matrix Uxi, 
and the sneutrino and charged slepton mixing matrices Uxi and Ux^. The chargino mixings 
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matrices are defined by the relation 

Mf'^s = {OR)M,{OLf, (28) 

where 

M.= f ° ^ V where X = ( „ (29) 

and Mc^°'^ is diagonal. The sneutrino mixing matrix • and the charged slepton mixing matrix 
■ are defined by the relations 

(^mj^)d^as ^U'mj^U} , (m?±)'^-^ = [/^m?[/|, (30) 

where the matrices m|_,_ and m? are given by the sum of the MSSM contribution and the respective 
Dirac leptogenesis contributions in (jS)) . The neutralino mixing matrix Un is defined by the relation 



(31) 



where 



/Ml -Mz sin 6*^ cos /3 sin (9^ sin /3 \ 

M2 Mz cos 9w cos P —Mz cos 6^ cos /3 

— sin 0w cos /? cos 9w cos /3 — /i 

\ Mz s'mOyj sin (3 —Mz cos 6w sin (3 —fx / 



(32) 



B Boltzmann Equations Including Second Order Processes 



In [21 , we derived the system of Boltzmann equations for Dirac leptogenesis up to processes first 
order in the heavy, decaying particles <I>i and and showed that the dynamics of these fields could 
be expressed using only two equations: one for the lepton number L^^ in the heavy field sector and 
one for the abundance Y^^^ of the scalar component of the supermultiplet. Here, we improve 
upon our previous calculations by including terms second order in (f>i. These terms only appear in 
the equation for the Y^^ abundance, which becomes^'' 






^^'^""VZ' ) ~^ ^ + 27ii ^-7a(^^-1), (33) 

where Ly^^ is the lepton number stored in sneutrinos, Lagg is the aggregate lepton number stored 
in the other, interacting fields that carry lepton number, and all abundances are normalized with 



^"Since the dynamics of scalar and fermion fields in $i and <&i are assumed to be the same, contributions of the 
form (j)\cj>i ij can be incorporated into 7a. 
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respect to the entropy density s. The reaction densities 71), 7^^, and 7]^ are given by 

'KiizY 



ID 



1 n 



eq 



sY2Vd 



Kijz) 
7 \K2{z) 

Kijz) 
K2{z) 



n 



eq 



(34) 
(35) 
(36) 



in terms of the Bessel functions Ki{z) and K2{z), the photon number density n^, and the rates 
Td, ^l, and Tr, which respectively represent the total decal width of the scalar component field 
01 and the individual decay widths for the processes (j) ^ £x and cj) — > i^^H^. The reaction density 
7A, which is associated with second order processes of the form ^i^i ij and (j)i4>i ij, is given 
by 

(37) 



647r4 



where T is temperature, s is the usual Mandelstam variable, and a{s) is the total reduced cross 
section for annihilations of </)i(/>i and (/>i(/>i into light fields. This is defined by the formula 



(38) 



where both t and s denote the Mandelstam variables. The limits of integration are given by 
t± = M"^^ — s{l =F '^)/2, with r defined below. 

In supersymmetric Dirac leptogenesis, the total reduced cross-section 7^, including all relevant 
decay processes, is 



- tot 
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(39) 



where x = s/M^^^ r = y^l — 4/x, and §2 and gy are the SU{2) and U{1)y coupling constants. 
The calculation can also be performed for the non-supersymmetric case, where the result is 
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(40) 



The effect of such second order annihilation processes on the parameter space of Dirac leptogenesis 
is shown in figure |3 where, for comparison, we show two sets of leptogenesis exclusion contours: 
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one representing no second-order processes and one representing annihilation in a super symmetric 
model. It is evident from this graph that second order processes do indeed lower the upper exclusion 
contour, though the effect is not a dramatic one. 
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Figure 4: This figure illustrates the effect of second-order processes of the form (picpi — > ij and 
<Pi(pi ij on the exclusion contours from leptogenesis. Contours are displayed for the case without 
annihilation (dashed line) and with annihilation (solid line). 
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